Thermal conductivity evidence for d x 2_ y 2 pairing symmetry 
in the heavy-fermion Celrln 5 superconductor 
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Quasi-two dimensional Celrlns contains two distinct domes with different heavy fermion super- 
conducting states in its phase diagram. Here we pinned down the superconducting gap structure of 
Celrlns in the second dome, located away from the antiferromagnetic (AF) quantum critical point, 
by the thermal transport measurements in magnetic fields H rotated relative to the crystal axes. 
Clear fourfold oscillation was observed when H is rotated within the afo-plane, while no oscillation 
was observed within the bc-plane. In sharp contrast to previous reports, our results are most consis- 
tent with d x 2_ y 2 symmetry, implying that two superconducting phases have the same gap symmetry 
which appears to be mediated by AF spin fluctuations. 

PACS numbers: 74.20.Rp,74.25.Dw,74.25.Fy,74.70.Tx 
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Superconductivity in heavy-Fermion (HF) compounds 
continues to be a central focus of investigations into 
strongly correlated electron systems. The relationship 
between the magnetism and unconventional supercon- 
ductivity, whereby the gap function A(fe) has nodes on 
the Fermi surface where A(fc)=0, is a particularly im- 
portant theme of research [lj. Many analyses have fo- 
cused on the scenario of superconductivity (SC) medi- 
ated by low energy magnetic fluctuations, often in prox- 
imity to a quantum critical point (QCP), where magnetic 
ordering temperature is driven to zero by an external pa- 
rameter such as pressure or chemical substitution. In- 
deed, unconventional SC appears in the vicinity of an 
antiferromagnetic (AF) QCP in most Ce-based HF com- 
pounds, including Celn 3 , CePd2Si2 [2(, CeCoIn 5 [3| and 
CeRhlns [4j , as well as organic and high-T c superconduc- 
tors. 

Notable counter examples have been recently reported 
in two Ce compounds, in which two distinct domes of 
different HF superconducting phases appear as a func- 
tion of pressure or chemical substitution. The first ex- 
ample is the prototypical CeCu2Si2, with one SC dome 
at low pressure around the AF QCP, and another dome 
emerging at high pressure distant from the QCP 0. The 
superconductivity in the low pressure dome is consistent 
with the magnetically mediated pairing. On the other 
hand, Cooper paring resulting from the Ce-valence fluc- 
tuations was proposed for the high pressure region, where 
no discernible AF fluctuations are found 

The second example is a quasi-two dimensional (2D) 
Celrlns (T c =0. 4 K) [3|, whose phase diagram is shown in 
the inset of Fig. I p], In this system the Rh sub- 
stitution for Ir increases the c/a ratio, acting as a neg- 
ative chemical pressure that increases AF correlations. 
In CeRhi_ a; Ir a ;In5, the ground state continuously evolves 
from AF metal (AFM) (x < 0.5) to superconductivity 
(x > 0.5). T c shows a maximum at x ~ 0.7 and exhibits 
a cusp-like minimum at x ^0.9, forming a first dome 



(SCI). The strong AF fluctuations associated with the 
AF QCP nearby are observed in SCI 0,0. In Celrln 5 
(x=l), T c increases with pressure and exhibits a max- 
imum (T c =l K) at P ~3 GPa, forming a second dome 
(SC2). The AF fluctuations are strongly suppressed with 
pressure, in SC2, far from the AF QCP @,El3- More- 
over, the nature of the AF fluctuations in magnetic fields 
in Celrlns is very different from that in CeCoIns and 
AF CeRhIri5 11, 12, T^. Thus the superconductivity in 
Celrlns at ambient pressure may be distinct from that 
in CeCoIns and CeRhlns, although all three compounds 
share similar quasi-2D band structure [14 1. 

Hence a major outstanding question is the nature of 
the microscopic pairing interaction responsible for the su- 
perconductivity in Celrlns. By analogy with CeCu2Si2, 
Celrln 5 in SC2 was suggested to be a strong candidate for 
the Ce-valence fluctuation mediated superconductor [||. 
To elucidate the pairing interaction, the identification of 
the superconducting gap structure is of primary impor- 
tance. Measurements of nuclear quadrupole resonance 
relaxation rate_E [l(| , thermal conductivity [H, [l6| and 
heat capacity [15| revealed that the superconductivity of 
Celrlns is unconventional, with line nodes in the gap. 
Very recently, from the measurements of the anisotropy 
between the interplane and intraplane thermal conduc- 
tivity k, the gap function of Celrln 5 was suggested to be 
of hybrid ty pe, k z (k x + ik y ), or E g symmetry [l6| . simi- 
lar to UPt3 [17], and in sharp contrast to the d x 2_ y 2 gap 
in CeCoIn 5 (and most likely in CeRhlns under pressure) 
[lit . However, as pointed out in Ref. [HI, the anisotropy 
of k alone is not sufficient to establish the hybrid gap, 
and further experiments aimed at clarifying the shape of 
the superconducting gap are strongly required. 

In this Letter, to shed light on the pairing mechanism 
of Celrlns, we performed the thermal transport measure- 
ments in magnetic fields H rotated relative to the crystal 
axes. We provide strong evidence that the gap symmetry 
is d x 2_ y 2 of Big symmetry. These results put a constraint 
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FIG. 1: (Color online) Temperature dependence of k/T as a 
function of T 2 in zero field and above H c c2 {H=2 T). Solid 
line represents k/T — kqq/T + AT 2 ; dashed line shows k/T 
obtained from the resistivity using Wiedemann-Franz law. In- 
set: phase diagram of Celrlns as a function of Ir concentration 
and pressure, from Refe.@,@- 
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FIG. 2: (Color online) Angular variation of the thermal 
conductivity, Ak = k — Ko, at T=200 mK with rotating 
H (a) within the &c-plane as a function of polar angle 6 
( | JJ" |=0.05 T) and (b) within the 2D ab-plane as a function 
of azimuthal angle 4> (|iT[=0.10 T). Here «o is the angular 
average of k. The thermal current q is applied parallel to the 
a-axis. The dashed line shows the twofold (cos 2(j>) variation. 
The solid line is the fit by «(</>) = Ko + C24, cos 24> + Ca4> cos 4<j>, 
where KojC^, and are constants. 



on the pairing mechanism in Celrlns. 

Single crystals were grown by the self-flux method. 
The bulk transition temperature is 0.4 K, and upper criti- 
cal fields parallel to the afe-plane and the c-axis, and 
ff c c 2 , are 1.0 T and 0.5 T at T=0 K, respectively. We 
measured the thermal conductivity k along the tetrago- 
nal a-axis (heat current q\\ a) on the sample with a rect- 
angular shape (2.8 x 0.45 x 0.10 mm 3 ) by the standard 
steady state method. To apply H with high accuracy 
(misalignment of less than 0.05°) relative to the crystal 
axes, we used a system with two superconducting mag- 
nets generating H in two mutually orthogonal directions 
and dilution refrigerator equipped on a mechanical rotat- 
ing stage at the top of the Dewar. 

Figure 1 depicts the temperature dependence of k/T 
as a function of T 2 at H — and in the normal state 
above H^ 2 . The overall behavior of k/T is similar to 
those reported in Refs. [l5| and [l6j]. In zero field, k/T 
decreases with decreasing T after showing a broad max- 
imum at T c , similar to UPt 3 [H and CePt 3 Si |21] . The 
value of k/T at Tr is nearly 30 % smaller than that re- 
ported in Ref. [la ] and nearly three times larger than 
that reported in Ref. [HI]. The dashed line is k/T ob- 
tained from the Wiedemann-Franz law, k/T = La/p, 

where Lq = ^- (^f-) 2 is the Sommerfeld value and p is 
the measured normal state resistivity. The observed k/T 
is close to Lq / p, which indicates that the heat transport 
is dominated by the electronic contribution. 

We first discuss the thermal conductivity in zero field. 
As shown in Fig. 1, k/T at low T is well fitted by 
k/T = kqq/T + AT 2 . The presence of a residual term, 



kqq/T, as T —> K is clearly resolved. This term in- 
dicates the existence of a residual normal fluid, which 
is expected for nodal superconductors with impurities. 
For a gap with line nodes, the magnitude of kqo/T is in- 



dependent of impurity concentration [22j. The present 
value of Knp/T ~ 2.2 W/K 2 m is close to that reported 
in Ref. [lfj, and, as discussed there, is in good agree- 
ment with the theoretical estimate for a superconductor 
with line node. Thus our thermal conductivity results 
are also consistent with the presence of line node in the 
gap function. 

The next important question is the nodal topology. 
Thermal conductivity is a powerful directional probe of 
the nodal structure: Recent measurements of both k 
and the heat capacity with H applied at varying ori- 
entation relative to the crystal axes established the su- 
perconducting gap structure in fc-space in several sys- 
tems 3 23, 24 1 . In contrast to fully gapped supercon- 
ductors, the heat transport in nodal superconductors is 
dominated by delocalized near- nodal quasiparticles. Ap- 
plied field creates a circulating supercurrent flow v s (r) 
associated with vortices. The Doppler shift of the en- 
ergy of a quasiparticle with momentum p, [E(p)— > E(jp)- 
v s -p], is important near the nodes, where the local en- 
ergy gap is small, A(p) < |i> s -p|- Consequently the 
density of states (DOS) depends sensitively on the an- 
gle between H and nodal directions (25)]. Clear twofold 
or fourfold oscillations of thermal conductivity and heat 
capacity associated with the nodes have been observed in 
UPd 2 Al 3 [26|, YBa 2 Cu 3 07 [13, CeCoIn 5 [Hill and K ~ 
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(BEDT-TTF) 2 Cu(NCS) 2 [29J when H is rotated relative 
to the crystal axes. 

Figures 2(a) and (b) show the angular variation of 
the thermal conductivity at 200 mK (k B T/A ~ 0.2) 
as H is rotated within the 6c-plane at |ff|=0.05 T 
(H/H£ 2 (T) ~0.14) and within the 2D a6-plane at 
|if|=0.1 T (H/H$(T) -0.14), respectively. Here 
=(H ,c) and 4> = (H,a) are the polar and azimuthal 
angles, respectively. For the field rotated within the ab- 
plane (9 = 90°) k(</>) exhibits a distinct oscillation as a 
function of <f>, which is characterized by peaks at (f> = 0° 
and ±90° and minima at around ±45°. As shown by 
the solid line, k(4>) can be decomposed into three terms, 
n((f>) = k + K 2( j ) + K44,, where k 2 = C 2 0COs2</> and 
K40 = C40 cos 40 have the two and four fold symmetry 
with respect to 0, respectively. We note that, as shown 
by the dashed line, k(<P) with minima at ±45° and peaks 
at ±90° cannot be fitted only by K 2 0-term, indicating 
the presence of the fourfold term. In sharp contrast to 
H rotating within the afe-plane, no oscillation is observed 
when rotating H within the 6c-plane; the amplitude of 
the oscillation is less than 0.2% of K n if it exists, where 
n n is the normal state thermal conductivity measured 
above i? c2 . The k 2 ^ term arises from the difference be- 
tween transport parallel and perpendicular to the vor- 
tices. Since for H within the &c-plane the field is always 
normal to q, K24, term is absent for this geometry. 

We address the origin of the fourfold oscillation. Fig- 
ures 3(a)-(d) display normalized by K n at 200 mK af- 
ter the subtraction of kq and K 2< p below _ff° 2 b (T) (~ 0.7 T). 
In the normal state above H^, no discernible oscilla- 
tion was observed (not shown). At 0.69 T just below 
H^iT), exhibits a minimum at cf) = 0° (C^ < 0). 
At H=0.5 T, oscillation vanishes. Further decrease 
of H leads to the appearance of distinct oscilla- 
tion that exhibits maximum at <fi = 0° (C40 > 0) at 
H^O.l and 0.25 T. Figure 3(e) shows the //-dependence 
of Ci^/Kn. There are two possible origins for the four- 
fold oscillation: (i) the nodal structure and (ii) in-plane 
anisotropy of the Fermi surface and Hgj. It should 
be stressed that the sign of just below is the 
same as that expected from the in-plane anisotropy of 
H c2 {Hc2 || (100) > H c2 || (110)) [U, whereas its 
sign at low fields is opposite. This immediately indi- 
cates that the origin of the fourfold symmetry at low 
fields is not due to the anisotropy of the Fermi sur- 
face or H®2 ■ Rough estimate of the amplitude of the 
fourfold term in layered <i-wave superconductors yields 
C H /K n = 0.082 ^j;^ ]n(y/32A/TrhT), where A is the 
superconducting gap, T is the QP relaxation rate, vf and 
v' F are the in-plane and out-of-plane Fermi velocities 31 1 . 
Using T - 1.3 xlO 11 s-^A/fcsTc ~ 5,v F - lxlO 4 m/s, 
and v' F ~ 5 x 10 3 m/s [§, [28| gives C^jK n ~2%, of the 
same order as the data. These results lead us to conclude 
that that the fourfold symmetry at low fields originates 
from the nodal structure. 
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FIG. 3: (Color online) (a)-(d) The fourfold component, 
normalized by k„ with H rotated in the ab plane at 
T=200 mK at |if|=0.69, 0.5, 0.25, and 0.1 T, respectively. 
The upper critical field is ^0.7 T at this temperature. 
The solid lines represent the fit by cos 40. (e) C^/nn at 
T=200 mK as a function of H/H$. (f) C^/K n at H=0.1 
and 0.25 T as a function of T/T c . 



The distinct fourfold oscillation within the afo-plane, 
together with the absence of the oscillation within the 
fee-plane, definitely indicates the vertical line nodes per- 
pendicular to the afe-plane, and excludes a horizontal line 
of nodes at least in the dominant heavy electron bands. 
Recall that in UPd 2 Al3 with horizontal line node, clear 
oscillations of k(9) are observed when rotating H within 
the ac-plane [26j. One could argue that the absence of 
oscillations within the fee-plane shown in Fig. 2(a) is due 
to relatively high temperature (fceT/A ~0.2), but the si- 
multaneous observation of the fourfold oscillation within 
the afe-plane at the same temperature rules out such a 
possibility. 

Thus the superconducting symmetry of Celrlns is nar- 
rowed down to either d x i- y i or d xy . Further identi- 
fication relies on the evolution of the oscillations with 
temperature and field. In the low-T, \ow-H limit, the 
Dopplcr shifted DOS shows a maximum (minimum) 
when H is along the antinodal (nodal) directions. How- 
ever, according to recent microscopic calculations, the 
pattern is inverted at higher T, H due to vortex scat- 
tering, and the fourfold components of the specific heat 
and of the thermal conductivity have similar behavior 
across the phase diagram 32J, [33|. In Fig. 3(f) we plot 
CV/ftn as a function of temperature. At iJ=0.1 T, the 
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FIG. 4: (Color online) H — T phase diagram for the fourfold 
term, H4<p, in H rotated within the 2D afo-plane. The solid 
(open) circles indicate the points where K40 shows a maximum 
(minimum) at <f> — 0°. The shaded regions correspond to the 
region where K40 has a minimum at the nodal direction, while 
g| has a maximum at the nodal directions outside this region 



netic propagation vector Q—(j, |, Q z ), has different Q z 
in these two phases. In fact, a remarkable difference be- 
tween CeCoIns and Celrlns is that the incommensurate 
AF order with Q^=0.298 strongly suppresses the super- 
conductivity in CeRhi-zCoaJns [351 ] . while they coexist 
in CeRhi_ x Ir a; In5 (3r3 |. 

In conclusion, the measurements of the thermal con- 
ductivity under rotated magnetic fields provide a strong 
evidence that the superconducting gap of Celrlns at am- 
bient pressure has vertical line nodes, and is of d x 2_ y 2 
symmetry. This indicates that two distinct domes of HF 
superconducting phases possess the same superconduct- 
ing symmetry, in which AF fluctuations appear to play 
an important role. The determined gap symmetry in 
Celrlns remote from the AF QCP further restricts theo- 
ries of the pairing mechanism. 

We thank H. Ikeda, K. Kontani, P. Thalmeier, and 
S. Watanabe for helpful discussions. 



sign change indeed occurs at T/T c ~ 0.25. Figure 4 dis- 
plays H — T phase diagram for the fourfold component. 
The solid (open) circles represent the points at which ob- 
served exhibits a maximum (minimum) at <j) = 0°, 
and the shading indicates the calculated anisotropy of the 
thermal conductivity for a cZ-wave superconductor from 
Ref. [33J. The shaded (unshaded) regions correspond to 
the region where has a minimum (maximum) for the 
field in the nodal direction. The calculation was done for 
a corrugated cylindrical Fermi surface, similar to that of 
the main FS sheet of Celrlns , and the results well repro- 
duce the observed sign change of k^. Since the minimum 
(maximum) of occurs at (j) = 45° inside (outside) 
the shaded region, the nodes are located at ±45°. We 
thus pin down the gap symmetry of Celrlns as d x 2_ y 2. 
While reconciling the existence of vertical lines of nodes 
with the results of Ref. 16] requires deviations from the 
perfect cylindrical symmetry of the Fermi surface [l9| . 
in the absence of detailed calculations for the realistic 
band structure of Celrlns (which would be desirable), 
the agreement with the computed phase diagram is re- 
markably good. 

The d x 2_ y 2 symmetry implies that the superconduc- 
tivity is most likely to be mediated by the AF spin fluc- 
tuations, not by the Ce- valence fluctuations [34} . Our 
result is also at odds with the hybrid gap function with 
horizontal node proposed in Ref. (l6| . To our knowl- 
edge, Celrlns is the first Ce-based HF compound in which 
d x 2_ y 2 symmetry is realized in a distinct superconducting 
phase remote from the AF QCP. It is intriguing that the 
superconductivity in SCI and SC2 phases has the same 
gap symmetry. A possible explanation for the double 
dome structure is that the active area on the Fermi sur- 
face for the superconductivity, which is nested by mag- 
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